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Abstract 

We present a measurement of the K-tt form-factor parameters for the decay Kl — > vr^i/. 
We use 328 pb~^ of data collected in 2001 and 2002, corresponding to ~ 1.8 million 
Kii'i events. Measurements of semileptonic form factors provide information about 
the dynamics of the strong interaction and are necessary for evaluation of the phase- 
space integral /j^ needed to measure the CKM matrix element \Vus\ for Kl vr^z^ 
decays and to test lepton universality in kaon decays. Using a new parameterization 
for the vector and scalar form factors we find A+=(25.6 ± 0.4stat. i 0.3syst.) x 10"'^ 
and Ao=(14.3 ± l.Tgtat. =t l.lstat.) x 10"'^. In the more usual quadratic expansion of 
the form factor the above result is corresponds to X'^ = A+, A'_|. = 2A^, Ag = Aq and 
Aq = (Aq + 0.000416)/2. Our results, together with recent lattice calculations of /jr, 
fx and /(O), satisfy the Callan-Trieman relation. 



1 Introduction 

Semileptonic kaon decays, Kl —>■ TT^i^iy, (Fig. 1) offer possibly the cleanest way 
to obtain an accurate value of the Cabibbo angle, or better, Vus- Since K —>■ n is 
a 0^ ^ 0^ transition, only the vector part of the hadronic weak current has a non 
vanishing contribution. Vector transitions are protected by the Ademollo-Gatto 
theorem against SU(3) breaking corrections to lowest order in (or rUg — Triu^d)- 
At present, the largest uncertainty in calculating Vus from the decay rate, is due 
to the difficulties in computing the matrix element {'n'\J^^'^\K). Using the notation 
of Fig. 1, Lorentz invariance requires that the above matrix element have the form 

< vrl >= {{P + PhUit) + {P- pUAt)) (1) 
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Fig. 1. Amplitude for Kl 
structure. 

where t = [P—pY 



■K^t^v. The gray region indicates the K — > t:W vertex 
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■2MEt^ is the only L-invariant variable. 



The form factors, FF, and f-{t) account for the non point like structure 

of the hadrons and the values of the FF's at t = differs from one because 
of SU{'i) corrections, i.e. because pions and kaons have different structure. The 
P — p = k + k' term in the /_ form factor, acting on the lepton term gives the 
lepton mass and is therefore negligible for K^s decays. The /„ form factor must 
be retained for i^^s decays. It has become costumary to introduce a scalar form 
factor fo(t) according to 



{nip)\u^^s\KiP)) = /(O) (P + p)^ Uit) + (P - p)^ /o(t) 



A 



hit) 



A 



with Aktt = M"^ — m?. The /+ and /o FFs must have the same value at t = 0. We 
have therefore factored out a term /(O). The functions and /o(t) are both 

unity at t = 0. If the FF are expanded in powers of t up to as 



t 



/+,o(t) = 1 + A'+,o -, + i,K 



,0 [ 2 



(2) 



four parameters: A'_,_, A'|, A'q and Aq need to be determined from the decay spec- 
trum in order to be able to compute the phase space integral which appears in 
the formula for the partial decay width. The problems with the four parameters 
above is the large correlations, in particular —99.96% between Ag and Ag and 
—97.6% between A'_,_ and A'^. It is not therefore possible to obtain meaningful 
results for the scalar FF parameters, see the appendix, eq. A.l. 

It is experimentally well established in K^^^ decays [1-3], that the vector form 
factor is equally described by a pole form: 



which expands to 1 + t/My + (t/MyY, neglecting power of t greater than 2. 

Recent results on Ke3 show that the vector form factor is dominated by the closest 
vector (qq) state with one strange and one light quark (or K-tt resonance in an 
older language) and are in good agreement with the results from fitting with 
a vector from factor as in eq. 2 [4]. There is however better consistency 

between the pole than the quadratic expansion fits, due mostly by the additional 
fluctuation introduced by the correlation, —95% between A', and A'f . The results 
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are also consistent with predictions from a dispersive approach [5,6]. We will 
therefore mostly use the following form for the vector form factor: 



Uit) = l + K—^+\l(-^) . (4) 

Kfj^3 decay pion spectrum measurements, reported in [1,7,8], have no sensitivity to 
Aq, see apendix. Therefore, all authors have fitted for a linear scalar form factor: 

fo{t) = 1 + Ao^. (5) 

Because of correlation this leads to incorrect answers for the value of Aq which 
comes out of the fit increased by ~3.5xthe coefficient of the term. To clarify 
this situation it is necessary to obtain a form for fo{t) with t and terms but 
with only one parameter. The Callan-Treiman relation [9] fixes the value of scalar 
form factor at t = A^n (the so called Callan-Trieman point) to the ratio of the 
pseudoscalar decay constants fx/ f-K- This relation is slightly modified by SU(2) 
breaking corrections [10]: 

/o(Ax.) = YJJo^ + ^CT, Act ^ -3.4 x 10"=^ (6) 

A recent parametrization for the scalar form factor [11] allow to take into account 
the constraint given by the Callan-Treiman relation: 

f,{t) = exp(-^\ogC-G{t)) (7) 



A 



where G{t) is obtained using a dispersion relation subtracted at t = Ak-k, such 
that C = /o(Ax7r)- As suggested in [11], a good approximation to eq. 7 is 

~ , , ^ , t A2 + 0.000416 / t \2 

This result is quite similar to ref. [12]. With KLOE, we can measure the pion 
energy spectrum {t = + — 2MEtt) spectrum since the value of Ki mo- 
mentum is known at a 0-factory. vr — /i separation is however very difficult at 
low energy. Attempts to distinguish pions and muons result in a loss of events of 
more than a factor of 2 and introduce severe systematic uncertainties. Therefore 
we use the neutrino spectrum that can be obtained without vr — yU identification. 

2 The KLOE detector 

The KLOE detector consists of a large cylindrical drift chamber (DC), sur- 
rounded by a lead scintillating- fiber electromagnetic calorimeter (EMC). A su- 
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perconducting coil around the calorimeter provides a 0.52 T field. The drift 
chamber [13] is 4 m in diameter and 3.3 m long. The momentum resolution is 
c"p^/p± ~ 0.4%. Two-track vertices are reconstructed with a spatial resolution of 
~ 3 mm. The calorimeter [14] is divided into a barrel and two endcaps. It covers 
98% of the solid angle. Cells close in time and space are grouped into calorime- 
ter clusters. The energy and time resolutions are ge/E = 5.7%/y^~(GeV) 

and cTr = 57 ps/^jE (GeV) © 100 ps, respectively. The KLOE trigger [15] uses 
calorimeter and chamber information. For this analysis, only the calorimeter sig- 
nals are used. Two energy deposits above threshold {E > 50 MeV for the barrel 
and E > 150 MeV for the endcaps) are required. Recognition and rejection of 
cosmic-ray events is also performed at the trigger level. Events with two energy 
deposits above a 30 MeV threshold in the outermost calorimeter plane are re- 
jected. 

3 Analysis 

Candidate events are tagged by the presence of a Ks tt+tt" decay. The 
tagging algorithm is fully described in Refs. [16] and [17]. The Kl momentum, 
PKl, is obtained from the kinematics of the — > KsKl decay, using the recon- 
structed Ks direction and the known value of p^. The resolution is dominated 
by the beam-energy spread, and amounts to about 0.8 MeV/c. The position of 
the production point, x,^, is determined as the point of closest approach of the 
Ks path to the beam line. The Ki line of flight {tagging line) is given by the K^ 
momentum, px^ = P0 — PKs ^"^^ the production vertex position, x^. All rele- 
vant tracks in the chamber, after removal of those from the Ks decay and their 
descendants, are extrapolated to their points of closest approach to the tagging 
line. 

For each track candidate, we evaluate the point of closest approach to the tagging 
line, Xc, and the distance of closest approach, dc- The momentum pc of the track 
at Xc and the extrapolation length, l^, are also computed. Tracks satisfying dc < 
avxy + b, with a = 0.03 and 6 = 3 cm, and —20 < /c < 25 cm are accepted 
as Kl decay products, r^y the distance of the vertex from the beam line. For 
each charge sign we chose the track with the smallest value of dc as a Kl decay 
product and from them we reconstruct the decay vertex. The combined tracking 
and vertexing efficiency for is about 54%. This above value is determined 
from data as described in Ref. [16]. Events are retained if the vertex is in the 
fiducial volume 35 < r^y < 150 cm and \z\ < 120 cm. 

Background from Kl^t^^t^^ , tt^tt^tt^ is easily removed by loose kinematic cuts. 
The largest background is due to Kl — > tt^'^c^z/ decays, possibly followed by early 
TT — > /ie decay in flight. For all candidate K^s events we compute A^e, the lesser 
between \Emiss — Pmissl assuming the decay particles are vre or fie. We retain 
events only if this variable is greater than 10 MeV. After the above kinematic 
cuts the efficiency for the signal is about 96% and the purity is about 80%. 
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A further cut on the distribution _E'jniss(7r''', /i )--E'miss(7r , fJ'~^) shown in Fig. 3 for 
Kl TTfiu and background events respectively, is apphed After the kinematic 



Fig. 2. £'jniss(7r^, ) versus £'miss(7r distribution from Monte Carlo. -^TTfiv 

(gray scale) and background (black dot). 

cuts described above, the contamination, dominated by Ki^nev decay is ~4%. 
Particle identification (PID) based on calorimeter information further reduces the 
contamination by ~2. 

Tracks are required to be associated with EMC clusters. We define two variables: 
diQ, the distance from the extrapolated track entry point in the calorimeter to the 
cluster centroid and (ix^TC; the component of this distance in the plane orthogonal 
to the momentum of the track at the entry position. We accept tracks with 
c^±,TC < 30 cm. The cluster efficiency is obtained from the KLOE Monte Carlo 
(MC), corrected with the ratio of data and MC efficiencies obtained from control 
samples. These samples, of 86% and 99.5% purity, are obtained from K^^, and 
i^es selected by means of kinematics and independent calorimeter information. 
The cluster efficiency correction versus Ey is shown in Fig. 3. 

For each decay track with an associated cluster we define the variable: Atj = 
^ci ~ ^ii {i = TT, e) in which td is the cluster time and is the expected time 
of flight, evaluated according to a well deflned mass hypothesis. We evaluate ti 
including also the time from the entry point to the cluster centroid [18]. We 
determine the e'^e~ collision time, to, using the clusters from the Ks- 

For this pupose, for each Kl decay track with an associated cluster, we deflne 
the variable: Atj = t^i —t^, {i = n, e) in which t^i is the cluster time and ti is the 
expected time of flight, evaluated according to a well-deflned mass hypothesis. 
The evaluating of tj includes the propagation time from the entry point to the 
cluster centroid [18]. 

An effective way to select the correct mass assignment, ire or en, is obtained by 
choosing the lesser of | At^+ — Ate- 1 and | At^- —Ate+\. After the mass assignment 




30 



6 



1 



0.99 
0.98 
0.97 
0.96 
0.95 
0.94 
0.93 







E,(MeV) -- 




20 40 60 80 100 120 140 160 180 200 



Fig. 3. Cluster efficiency correction versus E,^. 
has been made, we consider the variable 

-TLTOF — 



A'' - A' 
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Additional informations are provided by the energy deposition in the calorimeter 
and the cluster centroid depth. These measurements have been input to a Neural 
Network (NN). The value of -Rtof and that of the maximum of the NN outputs 
(NNmax) for the two charge hypotehsis are shown in Fig. 3 for data and MC. We 
retain events with i?TOF < l/6NNmax + 0.4 as indicated in Fig. 3. The resulting 




Fig. 4. NNmax vs i?TOF distribution (see text for the definitions) from Monte Carlo. 
Kl — >7r/iz^ (gray scale) and background (black dot). 



purity of the sample is 
IQQMeV (Fig. 3). 



97.5%, almost uniform in the fit range 21MeV < < 
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Fig. 5. Purity versus E,y. 

The form-factor parameters are obtained by fitting the E^, distribution of the 
selected events in the range 21 < E,y < 166 MeV. After subtracting the residual 
background as estimated from MC, we perform the fit using the following formula: 

dN 2° 

— (z) = NoY.A{z,j) X p(j, a;, A';, Ao) x e^otU) x FpsnU) (9) 
at .^-^ 

where p(j, A'^, A'_|., Aq) is the three-body differential decay width, and A{i,j) is the 
probability that an event with true value of E,y in the j^^ bin has a reconstructed 
value in the i^^ bin. The chosen bin size is 5.18 MeV, which corresponds to about 
1.7 (Jeu, where is the resolution on the neutrino energy. 

EpsR is the correction due to final state radiation. It is evaluated using the KLOE 
MC simulation, GEANFI [19], where FSR processes are simulated according the 
procedures described in Ref . [20] . FSR affects t-distribution mainly for high energy 
pions, i.e. for low t, where the correction is 3-5%. The slopes A'^, A'_[, and Aoare 
free parameters in the fit while the Nq constant is the total number of signal 
events. 



4 Systematic uncertainties 

The systematic errors due to the evaluation of corrections, data-MC inconsisten- 
cies, result stability, momentum mis-calibration, and background contamination 
are summarized in Tab. 1. 

We evaluate the systematic uncertainties on the tracking efficiency corrections by 
checking stability of the result when the selection of tracks is modified. We proved 
the validity of the method by comparing the efficiencies from data and MC control 
samples, and from the MC truth [18]. The uncertainty on the tracking efficiency 
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Source 


6X'^ X 10^ 


5A'| X 10^ 


6Xo X 10^ 


Tracking 


1.60 


0.47 


0.86 


Clustering 


1.50 


0.32 


1.45 


TOF + NN 


2.23 


1.16 


1.45 


p-scale 


1.10 


0.71 


0.81 


p-resolution 


0.61 


0.21 


0.01 


Total 


3.37 


1.50 


2.37 



Table 1 

Summary of systematic uncertainties on A^, X'^, Aq. 

correction is dominated by sample statistics and by the variation of the results 
observed using different criteria to identify tracks from Ki decays. Its statistical 
error is taken into account in the fit. We study the effect of differences in the 
resolution with which the variable dc is reconstructed in data and in MC, and 
the possible bias introduced in the selection of the control sample, by varying the 
values of the cuts made on this variable when associating tracks to Kl vertexes. 
For each variation, corresponding to a maximal change of the tracking efficiency of 
about ±10%, we evaluate the complete tracking-efficiency correction and measure 
the slope parameters. We observe changes of 1.6x10"^, 0.47x10"^, and 0.86x10"^ 
for A^ A" and Aq , respectively. We find a smaller uncertainty by comparing the 
efficiencies from data and MC control sample, and MC truth. However, we assume 
conservatively the changes in the result observed by varying the cut on dc as a 
systematic uncertainty. 

As for tracking, we evaluate the systematic uncertainties on the clustering effi- 
ciency corrections by checking stability of the result when the track-to-cluster 
association criteria are modified. Also in this case the uncertainty on the cluster- 
ing efficiency corrections is dominated by sample statistics and by the variation 
of the results observed using different criteria for the track-to-cluster association. 
We take into account its statistical error in the fit. The most effective variable 
in the definition of track-to-cluster association is the transverse distance, d±^TC- 
We vary the cut on (i_L,TC in a wide range from 15 cm to 100 cm, corresponding 
to a change in efficiency of about 19%. For each configuration, we obtain the 
complete track extrapolation and clustering efficiency correction and we use it to 
evaluate the slopes. We observe a corresponding changes of 1.5x10"^, 0.32x10""^ 
and 1.45x10"^ for A'_|_ A^ and Aq, respectively. 

We study the uncertainties on the efficiency of the and on the background evalua- 
tion by repeating the measurement on samples with modified PID and kinematic 
cut values, corresponding to a variation of the cut efficency from 90% to 95%. 
This allows to vary the background contamination from 1.5% to 4.5%. We observe 
a corresponding changes of 2.23x10"^, 1.16x10"^ and 1.45x10"^ for A'^ A'| and 
Ao, respectively. 
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The effect of the momentum scale and the momentum resolution have also been 
considered. We conservatively assume a momentum scale uncertainty of 0.1%, 
well above the known KLOE scale accuracy. We observe a corresponding changes 
of l.lxlO-^ 0.71x10-3 and 0.81 xlO'^ for A'^ A'| and Aq, respectively. 

We investigate the effect of momentum resolution by changing the value of the 
El, resolution by 3% as studied in [3]. The corresponding absolute changes are 
0.61x10-3, 0.21x10-3 and 0.01 xlO-^ for A' A^ amd Aq respectively. 



5 Results and interpretation 

About 1.7 Milion of K^^ events have been selected. We fit data, using the quadratic 
parametrization for the vector form factor and linear parametrization for the 
scalar form factor. The results are shown in Fig. 5. We obtain: 

a; = (22.3 ± 9.8stat. ± 3.4syst.) x lO^^ 
A'; = (4.8±4.9stat.±1.5syst.) X 10-3 
Ao = (9.1±5.W±2.4stat.) X 10-3 

with x^/dof = 19/29. The correlations are given by the matrix. 

Improved accuracy is obtained combining the above results with those from the 
Ke3 analysis [3]: 

AV = (25.5±1.5stat. il.Osyst.) X 10-3 
A'; = (1.4±0.7stat.±0.4sy,O X 10-3 



1 -0.97 0.81 
1 -0.91 
1 



We then find: 



A; = (25.6±1.5stat. ±0.9syst.) X 10' 

X 10"^ 
X 10' 



:i.5±0.7stat. ±0.4 



Ao = (15.4 ±1.8stat. ±1.1 



Stat. , 



^1 -0.95 0.29 ^ 



V 



1 -0.38 
1 



with x^/dof = 2.3/2 with the correlation given in the matrix on the right. 
Finally, we fit data using the parameterization for the scalar form factor given 
in eq. 8. Improved accuracy is obtained by using the pole parameterization for 
the vector form factor, truncated as in eq. 4. Dropping the " ' " indexes, we find 



A+ = (25.6 ± 0.4stat. ± 0.3syst.) x 10-3 

Ao = (14.3±1.7stat. ±l.lstat.) X 10 



,-3 
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Fig. 6. Residuals of the fit (top plot) and EnU distribution for data events superimposed 
on the fit result (bottom plot) 

with x^/dof = 2.56/3 and the correlation given at right. We remind the reader 
that terms are included as in eqs. 8 and 4. We note that using eq. 8, suggested 
in [11], the value of the phase space integral changes by only 0.04%. We find 
I{Kf,s) = 0.1026 ±0.0005. 

Finally, from the Callan-Treiman relation we compute /(O) = 0.964 ±0.023 using 
JkI f-K = 1.189 ± 0.007 from Ref. [21]. Our value for /(O) is in agreement with 
recent lattice calculations [22]. 

6 Conclusions 

A new measurement of the vr/iz/ form factors has been performed. Our result 

Kl^^ is in agreement with recent measurement from KTeV [1] and ISTRA± [8] 
and in disagreement with NA48 [2,7]. We also derive /(O) = 0.964 ± 0.023 in 
agreement with recent lattice calculations [22]. This agreement reinforces the 
credibility of the /(O) and Aq determinations. 
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A Error estimates 



It is quite easy to estimate the ideal error in the estimation of a set of parameters 
p=(pi, P2, ■ ■ ■ Pn) from fitting some distribution function to experimentally deter- 
mined spectrum. Let -F(p, x) be a probability density function, PDF, where p is 
some parameter vector, which we want to determine and x is a running variable, 
like t. The inverse of the covariance matrix for the maximum likelihood estimate 
of the parameters is given by [23]: 



92 InL 



dpidpj 

from which, for N events, it trivially follows: 



ij J F dpi dpj 

with dv the appropriate volume element. We use in the following the above re- 
lation to estimate the errors on the form factor parameters for one and two 
parameters expression of the form factors f+{t) and fo{t)- While the errors in any 
realistic experiment will be larger than our estimates, typically two to three times 
larger, it is still very important toward the understanding of the real problems in 
the determination of the parameters in question. 



A.l Kes decays 

For a quadratic FF, f{t) = 1 + X'^it/m"^) + {Xl/2){t/m^Y, the inverse of the 
covariance matrix G^^, the covariance matrix G+ and the correlation matrix 
are: 

5.937 13.867 \ 1 / 1.258^ -0.606 \ /I ■ 
13.867 36.2405 j ' N [ -0.606 0.509^ j ' \ 



N 




The square root of the diagonal elements of G+ gives the errors, which for one 
million events are 5A'^=0. 00126, 5A"=0. 00051. The correlation is very close to 
— 1, meaning a fit can trade A'^_ for A" and that the errors are enlarged. A fit for 
a linear FF, f{t) = 1 + X'_^{t/m'^) in fact gives A+=0.029 instead of 0.025 and an 
error smaller by ~3: 

5AV = ^(1g;^(1,1)) ' = 0.0004. 

A simple rule of thumb is that ignoring a term, increases A'^ by ~3.5xA'_|.. For 
decays the presence of a term in the FF is firmly established. It is however 
incorrect to try to fit for two terms connected by the simple relation A'^=2xA'^_2, 
both from theory and experiment. The above discussion justifies the use of eq. 4. 
The errors obtained above compare reasonably with the errors quoted in [1-3]. 
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A. 2 decays 



The scalar FF only contributes to K^^^ decays. Dealing with these decays is much 
harder because: a) - the branching ratio is smaller, resulting in reduced statistics, 
b) - the or t range in the decay is smaller, c) - it is in general harder to obtain an 
undistorted spectrum and d) - more parameters are necessary. This is quite well 
evidenced by the wide range of answers obtained by different experiments [1,7,8]. 
Assuming that both scalar and vector FF are given by quadratic polynomials as 
in eq. 2, ordering the parameters as Aq, A'J, \'_^ and A" , the matrices Gq]^^ and 



are: 



N 



/1.64 


5.44 


1.01 


3.90 \ 




( 


63.9^ 


-1200 


-923 


197 \ 


5.44 


18.2 


3.01 


12.3 


1 




-1200 


18.8^ 


272 


-59 


1.01 


3.01 


1.47 


4.24 


' N 




-923 


272 


14.8^ 


-49 


\3.90 


12.3 


4.24 


13.8/ 




\ 


197 


-59 


-48 


3.4V 



and the correlations, ignoring the obvious terms, are: 

/ -0.9996 -0.974 0.91 \ 

0.978 -0.919 . 
^ -0.976 ] 



(Al) 



All correlations are very close to —1. In particular the correlations between Aq 
and Aq is —99.96%, reflecting in vary large (5Aq and ^Aq errors. We might ask 
what the error on Aq and Aq might be if we had perfect knowledge of A^ and \"^. 
The inverse covariance matrix is give by the elements (1,1), (1,2), (2,1) and (2,2) 
of the Gq^_,_ matrix above. The covariance matrix therefore is : 

Go(A(„ A'q' for a;, A'I known) = 1 ^'^^ ^.42 ) • 

For one million events we have (5Ao=0.0024, about 4x the expected value of Aq. In 
other words Aq is likely to be never measurable. It is however a mistake to use a 
scalar form factor linear in t, because the coefficient of t will absorb the coefficient 
of a term, again multiplied by ~3.5. Thus a real value Aq=0.014 is shifted by 
the fit to 0.0161, having used eq. 8. Fitting the pion spectrum from 1 million 
i^^s decays for Aq, A'^ with the form factors of eq. 8 gives the errors (5Ag~0.001 
and 5A'_,_~0.0011. Combining with the result from a fit to 1 million K^,^ with the 
FF of eq. 4 for which 5A'+~0.0004 gives finally 5A'q~0.00082, 5A'+~0.00038 and a 
Aq-A'^ correlation of —29%. We hope to reach this accuracy with our entire data 
sample, ~5x the present one, and a better analysis which would allow using the 
pion spectrum. 
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